Using self-driven microswimmers for particle separation 



(N 

O 

(N 



in 



o 



o. 
>; 

On, 
On. 

o- 
in; 

o: 



X: 



W. Yangi'2, V. R. Misko\ K. Nelissen^^^^ M. Kong", and F. M. Peeters^'^ 
^ Departement Fysica, Universiteit Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium 
''College of Materials Science and Engineering, Taiyuan University of Science and Technology, Taiyuan 030024, P- R- China 

^ Departamento de Fisica, Universidade Federal do Ceard, 
Caixa Postal 6030, Campus do Pici, 60455-760 Fortaleza, Ceard, Brazil 
'^Institute of Plasma Physics, Chinese Academy of Sciences, Hefei 230031, P. R. China 

(Dated: March 16, 2012) 

Microscopic self-propelled swimmers capable of autonomous navigation through complex environ- 
ments provide appealing opportunities for localization, pick-up and delivery of micro-and nanoscopic 
objects. Inspired by motile cells and bacteria, man-made microswimmers have been fabricated, and 
their motion in patterned surroundings has been experimentally studied. We propose to use self- 
driven artificial microswimmers for separation of binary mixtures of colloids. We revealed different 
regimes of separation including one with a velocity inversion. Our finding could be of use for various 
biological and medical applications. 

PACS numbers: 36.40.Wa, 82.70.Dd, 87.15.hj 



I. INTRODUCTION 

Quasi-one-dimensional (QID) colloidal systems play 
an increasingly important role as model systems to study 
a variety of collective phenomena in low dimensional con- 
densed matter physics. The well-controlled way in real 
space and time, and the tunable interparticle interac- 
tion potential of colloidal systems lead to insights of a 
wide range of systems. Most of the studies on QID col- 
loidal systems focused on stripe structures and crystal 
defects , solid- liquid phase transitions [1, 0| , single- 
file diffusion (T HlOl an d transport behavior under exter- 
nal driving force |lll4l4l |. including, e.g., dragged par- 
ticles [11]. Transport properties of colloidal particles 
driven by an external force in narrow channels became 
a growing topic during the last decade. This is by part 
because the transport behavior might provide a deeper 
understanding of several biological systems, e.g., ion- 
channel transport in cell membranes [1^, , DNA ma- 
nipulations and separations 18-20j. Also it is a read- 
ily available physical model to mimic other types of 
particles, including different ions and electrons, which 
might lead to practical use in microdevices and nanode- 
vices [lllll^]. Moreover, the transport behaviors are im- 



portant in real systems of medicine, food production and 
field-response materials as well. There are also studies 
of transport properties of particles which arc driven in 
other shapes of patterned environments 23, 24]. 

In nature and technology, many systems are mixtures 
of different particle types. The distinction and interac- 
tion, mixture and separation of each particle type give 
rise to an exceedingly rich phenomenology as compared 
to monodisperse systems. This motivited a number of 
studies on complex mixtures and, in particular, binary 
mixtures (see, e.g., psl - ls^ )). 

One of the fundamental problems related to binary 
mixtures is their separation (see, e.g., (25| - |29| ). This 
problem is commonly tackled by using different asym- 
metric potentials for separating, e.g., interacting binary 



mixtures driven on periodic substrates [25|-|29j). or fer- 
rofluids [Slj. Mixtures of particles can be driven either 
along the symmetry axis of the potentials, or in the trans- 
verse direction [32-| which was demonstrated in the con- 
text of separation of macromolecules, DNA, or even cells 
(see, e.g., [H-ll^). 



Note that in the previous studies a driving force (e.g., 
gravity [1^ or external field [HI, [13, 113] ) applied to 
all the particles of the system (or to one of the species) , 
and the separation occurred either due to the different 
dynamic response of the species to the driving/potentials 
or due to their different self-diffusion constant. 



In contrast to that, here we propose a new mechanism 
of particle separation. It is based on injecting special 
particles which are able to move in the binary mixture 
and interact with its species. These can be either parti- 
cles driven by external forc es ( which describes a typical 
situation in microrheology |38|) or self-driven particles 
i.e., driven by various self-phoretic forces produced by 
a chemical (self-chemophoresis) (see, e.g., the recent re- 
view paper [S^), electrical (self-electrophoresis) or ther- 
mal (self-therniophoresis) gradient that the particle gen- 
erates around itself (see Ref . [40] and references therein) . 
For example, so-called Janus particles (/j,m-sized parti- 
cles covered by gold on one of the hemispheres) can move 
( "swim" ) when being illuminated by light (details of the 
mechanism of motion of these artificial "microswimmers" 
(MS) are beyond the scope of this work and can be found 
in [lo, |4l|). Due to their amphliphilicity, Janus particles 
are extensively studied as a model of structure-directing 
amphiphiles (see, e.g., [12]) and in relation to thei r ap - 
plications in new functional materials (see, e.g., [43ll44|). 
Here we demonstrate that artificial MS, when moving in 
a binary mixture, are capable of selectively driving and 
separating the flows of different species of the binary mix- 
ture. 
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II. SIMULATION 

We consider a binary system of paramagnetic colloidal 
particles confined in a 2D infinitely long narrow hard- 
wall channel in an external perpendicular magnetic field. 
The particles (including the microswimmer — for ex- 
perimental realization of paramagnetic microswimmers 
see Refs. [H, interact via a repulsive dipole-dipole 
potential (see, e.g., experiments [1, [11]) Vij{ri,fj) — 



Qi Qj 1 1 '^i 



where Qi = [XqI Att is the "effec- 



where the unit length ao = 10 /xm is on the order of 
the average distance between particles [1, [13] . Thus the 
inter-particle interaction can be characterized by the di- 
mensionless coupling parameter F = Q\/ a^ksT ^ where 
ks is the Boltzmann constant and T the ambient tem- 
perature. In our simulations, we choose T = 200 which 
corresponds to the strongly correlated regime, i.e., below 
the melting transition. Correspondingly, the unit force 
is Fq — TkBT/aQ. For the chosen F — 200 and room 
temperature, Fo = 8.28 • 10" N. 

The motion of colloidal particles is investigated us- 
ing the Langevin equations of motion in the overdamped 
regime, i.e., the Brownian dynamics (BD) method. This 
approach neglects hydrodynamic interactions as well as 
the short-time momentum relaxation of the particles (see, 
e.g., [H, IS, 53] ) • The overdamped equations of motion 
of particles are: 



— = ^{~Vr~yV,,(r, 



rj)+Fdr^ + F,{t)}, (1) 



where Di is the self-diffusion coefficient, which is cho- 
sen as Df — 1.0 fim? / s for type A (and for MS) and 
Df = 0.7 fim? / s for type B particles, according to the 
experiment [l^, EH . The first term of the rhs of Eq. ([T]) is 
the sum of all repulsive interaction forces acting on each 
particle, the second term is the "self-driving force" acting 
on/?/ on the MS (or "motor" particle L48,49i])- Recently, a 
similar approach has been used for the description of mo- 
tion of swimming bacteria and self-driven particles [50| . 
Note that force F^j.^ entering into the Langevin equa- 
tion ([T]) is an effective force which describes the propul- 
sion mechanism on average (see, e.g., the recent stud- 
ies [sil . Is^ on the Brownian motion of a self-propelled 
particle). In experiment, this effective "driving force" 
can be tuned by, e.g., changing the intensity of the inci- 
dent light (in case of Janus particles) , the geometry of the 
experiment [4ll |. or the number of MS particles, etc.; we 
do not discuss here the relation between F^rx and those 
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five charge". Mi is the magnetic moment, and ri is the 
coordinate of the «-th particle; /ig is the magnetic permit- 
tivity. We assume that there are two types of particles, 
A and B, which differ by size and charge. We call them 
"small" (A) and "big" (B) particles. The system con- 
sists of Na small particles with Qa, and Nb big parti- 
cles with Qb — 8Qa- The effective charge of self-driven 
MS is chosen Qd = Qa- We introduce a unit energy 

— Qa /^o °f tli*^ inter-particle interaction Vij (ri ,rj), Q 




FIG. 1: (a) Equilibrium configuration of a system with equal 
number of particles of the two species Na ~ 40, Nb ~ 40 and 
d ~ 8Qa- (b) The average velocity of particles of type A, 
(V^) (black open squares connected by solid lines) and B, 
{V^) (red dots connected by lines), as a function of Fdrx- 
The microswimmer (MS) is shown by a blue filled circle. The 
same symbols are used in Figs. 2 to 5. The inset shows a 
zoom of the inversion velocity region (2). 



parameters. Furthermore, without loss of generality we 
consider here only the x-componcnt of Fdr- This is justi- 
fied by the fact that the (average) direction of motion of 
a MS (i.e., a Janus particle illuminated by light) can be 
controlled by, e.g., direction of the incident light and by 
using different illumination patterns jiol . l4l| , or by addi- 
tionally applied drift force 41|. Furthermore, in our case 
the motion of a MS across the channel is restricted by 
the boundaries. For simplicity, we ignore rotational mo- 
tion of a MS [5^. The last term is the thermal random 
force that obeys the following conditions: {Fi{t)) — 0, 
and {F,^{t)F,pit')) = 2 {klT^/D,) 5{t ~ t')5,,5^0. 



As a simulation cell, we use a unit box of length 
Lx = 100 (Imm) and width Ly = 12 (120/Ltm) (the box 
dimensions are close to those in the experiments [slfl^). 
Periodic boundary conditions (PBC) are imposed in the 
X direction, and hard wall boundaries in the y direction. 
We use a cutoff (r = 6) for the dipole short-range inter- 
particle interaction. We also checked that by doubling 
the length of the simulation cell our results do not change. 

The initial equilibrium configurations of the system 
were calculated for Fdrx — 0. Typical equilibrium dis- 
tribution is shown in Fig. [Ija) for Na = Nb = 40 
and Nms = Ij which is consistent with our previous re- 
sults [30] (without MS particle; note that in our present 
simulations Nms ^ ^Ai Nb)- Then we calculated the 
average velocities of A and B particles as a function of 
Fdrx- Simulations were typically performed for 10^ time 
steps Ai = 10"^ (s). 
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III. SEPARATION OF THE PARTICLE FLOWS 

The average velocity (along the channel) (Vx) of parti- 
cles A and B versus Fdrx is shown in Fig. [ijb) for a typi- 
cal system with Na = Nb — 40. The main characteristic 
feature of the function {Vx){Fdrx) is a well-pronounced 
particle flow separation for a broad range of F^rx- Note 
that particles A (small) move faster than their counter- 
part. Depending on F^rx, we can distinguish four differ- 
ent regions (shown in Fig.[ljb) as (1) to (4)). 

Region (1): Rigid body motion 

This regime corresponds to very small Fdrx when the 
system is in a "rigid body" motion, i.e., velocities {V,f) 
and (Vj?) coincide. The system is in a solid state (F = 
200), and only elastic deformations are present [l3l |. 

Region (2): Inverse velocity motion 
This regime is observed for small F^rx when the veloc- 
ity of the MS is high enough to locally melt the solid, 
however, the velocites are inversed as compared to the 
broad region (3), i.e., particles A (small) move slower 
than particles B. This unusual behavior is explained by 
the stronger interaction of the MS with particles B than 
with particles A. As a consequence, particles B are car- 
ried along by the MS while the dynamical friction due 
to the particle motion is very small in this case and can 
be neglected. This regime is further illustrated in Fig. [2] 
The snapshots (Figs. UKa-c)) show that the system ad- 
justs itself to the slow "swimming" of the MS, without 
appreciable changes in the structure. However, the local 
density profile (Fig.[2ljd)) indicates that the local density 
of particles B in front of the moving MS is higher than 
that of particles A. This means that particles B (big) are 
carried along by the MS while particles A (which interact 
weaker with the MS) "slip" and thus move slower. 

Region (3): Strong flow separation 

This regime is observed for a broad range of intermedi- 
ate Fdrx'- this is the main and the most robust regime. 
With increasing velocity, the dynamical friction becomes 
increasingly important, and, as a result, type A particles 
(which are characterized by a larger self-diffusion coef- 
ficient) start to move faster. The function (y^){Fdrx) 
reaches its maximum at about Fdrx = 2 and then grad- 
ually decreases while the function {V^){Fdrx) increases 
only slowly in this region and then decreases. It is worth 
noting that in this regime, the function (yJI^){Fdrx) first 
rapidly decreases. This decrease is accompanied by a si- 
multaneous sharp increase in the function {y^){Fdrx)- 
The origin of this behavior is explained by the fact that 
the motion of the two species (i.e., particles A and B) be- 
comes less correlated as compared to that in regime (2). 
This loss of correlation in motion of the different species 
resembles the transition from elastic to plastic motion. 
The observed splitting of the average velocities is similar 
to the splitting in the velocities of adjacent layers found 
in a "vortex Wigner solid" [sj, [s^ where the transition 
to plastic motion was induced by shear stress. The snap- 
shots and the local density profiles shown in Fig. |3] illus- 
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FIG. 2; (a-c) Snapshots of the configurations of the binary 
system (small particles of type A are shown by empty black 
squares; large particles of type B are shown by red dots) in 
presence of a microswimmer (shown by a blue filled circle), 
for Na = 40, Nb = 40, Qd = 8Qa and Fdr^ = 0.4, which 
corresponds to Regime (2) (see Fig. 1) of inverse velocities. 
The snapshots are recorded at f = 3 (a), 60 (b), and 150 
(c). (d) The local density distribution of the particles versus 
distance to the MS, AX. 



trate this regime. In contrast to regime (2) (cp. Fig. [2]), 
the motion of the MS strongly influences the particle dis- 
tribution: in particular, we observe a pronounced effect 
of particle separation in real space (i.e., not only in ve- 
locity space), or the formation of clusters of particles A 
which "trap" the MS and follow it (see Figs, c)). 
The local density of particles A has a sharp peak in front 
of the MS (Fig. |3i;d)), due to the MS- induced compres- 
sion of the particle configuration in front of the MS. At 
the same time, the particles behind the MS are relaxed 
resulting in no peak in the local density. 

Region (4): Fast MS motion. 

This regime corresponds to large Fdrx ■ finally, the system 
undergoes a transition to a "quasi-rigid body" regime 
when the MS moves too fast through the binary solid to 
produce a responce in form of flow separation. 



In order to analyze the role of the width of the chan- 
nel, we performed simulations for wider channels, and we 
found that the effect of separation remains qualitatively 
similar: for a broad range of Fdr, the average velocity of 
small particles, {V^){Fdrx), is larger than that of large 
particles, {V^){Fdrx)- However, the more gentle effect of 
velocity inversion can disappear, depending on the con- 
centration of the particles. The role of imbalance in prop- 
erties of the particles is discussed in the next section. 




FIG. 3: (a-c) Snapshots of the configurations of the binary 
system in presence of a microswimmer (the same symbols are 
used as in Fig. 2), for Na = 40, Nb = 40, Qd = 8Qa and 
Fdrx = 1.2, which corresponds to Regime (3) (see Fig. 1) of 
strong flow separation. The snapshots are recorded at t — 3 
(a), 600 (b), and 23100 (c). (d) The local density distribution 
of the particles versus distance to the MS, AX. 



FIG. 4: (a-c) The average velocity (14) of particles A and 
B as function of Fdrx for a system with Na ~ Nb ~ 40 
and different charge of the MS, Qa = 2Qa (a), 12Qa (b), 
and 128(5 A (c). The inset in (a) shows the velocity of the 
driven aprticle versus Fdrx for systems with Qd ~ 2Qa and 
Qd ~ 128Qa- (d) The maximum ratio of the velocities of 
particles A and B, VRmax, as a function of charge of the MS, 
Qd. 



IV. ROLE OF IMBALANCE IN PARTICLE 
PROPERTIES 

Above, we revealed the general features of the sepa- 
ration of the particle flows by the MS "swimming" in a 
binary mixture and explained the mechanisms of the ob- 
served effect for the different separation regimes. Here we 
analyze the influence of the imbalance in particle prop- 
erties, i.e., differences in charge of the MS and the con- 
stituent particles and the relative density of particles, on 
the observed effect. 

We systematically examined systems with the charge 
of the MS varying in a broad range from Qd — Qa to 
Qd — 128Qa- The results of our calculation of (Vx) ver- 
sus Fdrx for Qd — 2Qa, 12(5a, and 128Qa are presented 
in Figs. Ufa-c). These results show that the position of 
the maximum average velocity (V^) shifts towards larger 
Fdrx with increasing the Qd/QA ratio. We explain this 
by the fact that increasing the interparticle interaction 
enables the system response to the motion of the MS 
even for large velocities, with simultaneous extension of 
the initial rigid-body region (see Fig. HJc)). Note that 
the region of inverted velocities exists only for a certain 
range of values of Qd/QA ^ 10. 

We analyzed the maximum effect of flow separation, 
i.e., the ratio {V^)/{V^) = VRmax and calculated its 
dependence on Qd/QA (Fig. [Hd)): the initial fast grow 
is followed by a wide "plateau" (for Qd/QA = 10 to 60) 
where {V^) / (yj^) ~ 4, then followed by a considerable 
increase for Qd/QA > 60. 

We also studied the role of the imbalance in the den- 
sity of particles A and B on the flow separation. We 
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FIG. 5: (a-c) The function (14) for particles A and B versus 
the driving force for a system with Qd ~ 8Qa and different 
relative densities of particles A and B. (d) The maximum ratio 
of the velocities of particles A and B, VRmax, as a function 
of number of particles A, A''^. 



found that the position of the maximum of (V^) shifts 
towards larger Fdrx with increasing the Na/Nb ratio. At 
the same time, the maximum effect of flow separation, 
(y^) / {Vj^) turns out to be a non-monotonic function of 
Na/Nb, with a maximum at Na/Nb ~ 2 (see Fig. [5]). 
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V. CONCLUSIONS 

We demonstrated that particles in a binary mixture 
can be effectively separated by using self-driven parti- 
cles, or microswimmers, that move through the system. 
We analyzed the main features of the particle separation 
and explained mechanisms of different regimes including 
one with a velocity inversion. Our theoretical findings 
can be verified in experiments with microswimmers in 
colloidal binary mixtures. Indeed, the revealed regimes 
of particle flow separation are based on a rather gen- 
eral interplay between the interaction and the dynamical 
friction. In addition, in our model we used realistic ex- 
perimental parameters. For example, our estimates for 
the "driving force" needed to observe the predicted be- 
havior is on the order of (0.1-0.2)pN which agrees with 
the estimates for light-driven Janus particles in the ex- 
periments (iol. l4l|. Furthermore, the same experimental 
studies reported on: (i) the high control of the direction 
of motion of microswimmers (i.e., by employing various 
illumination patterns) and (ii) the tunability of the "driv- 



ing force" of light-driven Janus particles (i.e., by tuning 
the intensity of the illumination). These two ingredients 
are of key importance for the possibility of experimental 
realization of the present theoretical findings. Moreover, 
the approach used in this work is also applicable to mi- 
crorheology. The results of our study can be useful for 
mixture separation in biology and medicine. 
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